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Summary 

This paper shows how a basic property of unitary transformations can be used for meaningful computations. This approach 
immediately leads to search-type applications, where it improves the number of steps by a square-root - a simple minded search 

that takes N steps, can be improved to approximately 4N steps. The quantum search algorithm is one of several immediate con- 
sequences of this framework. Several novel search-related applications are presented. 

1. Introduction Several interesting problems in computer science can be looked upon as search problems. There 
are two categories of such problems. First, where the search depends on data in memory - this is the database search 
kind of problems. Alternatively, the search could be based on a function known in advance - many NP-complete prob- 
lems and cryptography problems can be expressed in this form. For example the SAT problem of NP-completeness 
asks whether there exists a combination of binary variables that satisfies a specified set of Boolean equations - this 
can be looked upon as a search of the state space of the binary variables. In cryptography, the well-known 56-bit DES 

code (Data Encryption Standard) can be cracked by an exhaustive search of 2^ items [BBHT] [Phone]. 

It aroused considerable interest when it was shown that it was possible to improve upon the obvious classical 
bound for exhaustive search by resorting to quantum mechanics [Search] [BBHT] - the intuitive reason for this 
improvement was that quantum mechanical systems can be in multiple states and simultaneously explore different 
regions of configuration space. This improved the number of steps by a square-root, i.e. a simple minded search that 

takes N steps, could be improved to approximately Jn steps. The quantum search algorithm was derived using the 
Walsh Hadamard (W-H) transform and it appeared to be a consequence of the special properties of this transform. 
Subsequently [Gensrch] showed that similar results are obtained by substituting any unitary transformation in place 
of the W-H transform. This means that a variety of unitary transformations could be used in place of the W-H trans- 
form and this leads to algorithms for several different problems. This paper describes the approach of [Gensrch] and 
shows how it can be extended to solve various structured problems. 

2. Framework A function f(x), x = 0, 1, . . . (N - 1 ) , is given which is known to be zero for all x except the sin- 
gle point x = t , the goal is to find t (t for target). The obvious classical technique of searching by looking at the N 
values of x, one by one, would clearly take O(N) steps. 

Assume that we have at our disposal a unitary transformation U that acts on a system with N basis states. First 
map each value of x to a basis state and start with the system in the basis state s (s for start). If we apply U to s, the 
amplitude of reaching t is U , and if we were to make a measurement that projects the system into a unique basis 



2 

state, the probability of getting the right basis state would be t/ . It would, therefore, take Q. 



^ U tf J 



repetitions 



of this experiment before a single success. This section shows how it is possible to reach state t in only O 



1 



steps. This leads to a sizable improvement in the number of steps if | U (s \ « 1 . 

Denote the unitary operation that inverts the amplitude in a single state x by / . In matrix notation this is the 

diagonal matrix with all diagonal terms equal to 1 , except the xx term which equals -1 . 

th 

v denotes the column vector which has all terms zero, except for the x term which is unity. 

Consider the following unitary operator: Q = -I U I ( U, since U is unitary U 1 is equal to the adjoint, i.e. the 
complex conjugate of the transpose of U. We first show that Q preserves the two dimensional vector space spanned 

by the two vectors: v and (U *v f ) (note that in the situation of interest, when {/ is small, these two vectors are 
almost orthogonal). 

-1 T 

First consider Qv . By the definition of Q , this is: -I U LUv . Note that v v is an N X N square matrix 

A SIS JC A, 

T T 

all of whose terms are zero, except the xx term which is 1 . Therefore 1=1- 2v v. & 1=1- 2v v , it follows: 

(1) Qv s = -{I-2v s v T s )U~ l {I-2v/ t )Uv s = -{I-2v/ s )U~ X Uv s + 2{I-2v s v T s )U~ l {v/ t )Uv s 

-1 T 
Using the facts: U U = I and v v = 1 , it follows that: 

(2) Qv s = v s + 2(I -2v s v T s )U~\v t vJ)Uv s . 

T 

Simplify the second term of (2) by the following identities: v f Uv s = U & since U is unitary, its inverse is equal to 

T -1 * 

its adjoint (the complex conjugate of the transpose) v U v f = U f . 

(3) Qv s = v s (l-4\U ts \ 2 ) + 2U ts (U- 1 v t ) 

Next consider the action of the operator Q on the vector U ^v f . Using the definition of Q (i.e. 
Q = -I S U ^ IjU ) and carrying out the algebra as in the computation of Qv ? above, this yields: 

(4) Q(U-\) = -I s U- l I t U(U-\) = -I s U- l I t v t = I s U-\. 

T T — 1 * 

Writing / as / - 2v„v c and as in (3), v„ U v. = U t „ : 



(5) Q(U \) = (U l v t )-2v s v T s (U \) = (U \)-2[/* 4 , v 

It follows that Q transforms any superposition of v & (U *v ) into another superposition of the two vectors, 



-1 



thus preserving the two dimensional vector space spanned by v & ( U v ( ) . (3) & (5) may be written as 



(6) Q 
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U v. 



It follows as in [BBHT], that if we start with v , then after r) repetitions of Q we get the superposition 
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av+a.(U v.) where a s cos(2ri [/. J) & \a\ = sin(2ri \U f I) if \U. I « 1 . If ri = —r— — - , then we get the 
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superposition U 1 v { ; from this with a single application of U we can get v f . Therefore in o[ r-^— : I steps, we can 



ts\ 



start with the s-state and reach the target state t with certainty. 



3. Quantum Operations The interesting feature of the analysis of section 2 is that U can be any unitary trans- 
formation, whatsoever. Clearly, it can be used to design algorithms where U is a transformation on the qubits in a 
quantum computer - the object of this paper is to present some such applications. Quantum mechanical operations 
that can be carried out in a controlled way are unitary operations that act on a small number of qubits in each step. It 
is possible to design a variety of quantum mechanical algorithms using just a few elementary quantum mechanical 
operations. Two of the elementary unitary operations needed are: the W-H transformation operation and the selective 
inversion of the amplitudes of certain states. 

A basic single bit operation in quantum computing is the operation M - this is represented by the following 



matrix: M = — 
72 



1 1 
1 -1 



, i.e. a bit in the state is transformed into a superposition: ( — , — . Similarly a bit in state 

Li a/ -t-^ 



1 1 



1 is transformed into \-^=, — \=\ . In a system in which the states are described by n bits (it has N = 2" possible 

VV2 J2.J 

states) we can perform the operation M on each bit independently in sequence thus changing the state of the system. 

n n 

The state transition matrix representing this operation will be of dimension 2x2. Consider a case when the start- 
ing state is one of the 2 n basis states, i.e. a state described by an n-bit binary string with some Os and some Is. The 
result of performing the operation M on each bit will be a superposition of states described by all possible n-bit binary 



strings with the amplitude of each state being ±2 . To deduce the sign, observe that from the definition of the matrix 



M, i.e. M = — 

72 



1 1 
1 -1 



, that the phase of the resulting configuration is changed only when a bit that was previously a 



1 remains a 1 after the transformation is performed. Hence if x be the n-bit binary string describing the starting state 

and y the n-bit binary string describing the resulting string, the sign of the amplitude of y is determined by the parity 

_ x ' y 

of the bitwise dot product of x and y , i.e. (-1 ) . This transformation is the W-H transformation [DJ]. This opera- 
tion (or a closely related operation called the Fourier Transformation [Factor]) is one of the things that makes quan- 
tum mechanical algorithms more powerful than classical algorithms and forms the basis for most significant quantum 
mechanical algorithms. 

The other transformation we will need is the selective inversion of the phase of the amplitude in certain states. 
Unlike the W-H transformation and other state transition matrices, the probability in each state stays the same since 
the square of the absolute value of the amplitude in each state stays the same. The following is a realization based on 
[BBHT]. Assume that there is a binary function f(x) that is either or 1 . Given a superposition over states x, it is 
possible to design a quantum circuit that will selectively invert the amplitudes in all states where f(x) = 1 . This is 
achieved by appending an ancilla bit, b and considering the quantum circuit that transforms a state \x, b) into 
\x, f(x)XOR b) (such a circuit exists since, as proved in [Revers], it is possible to design a quantum mechanical cir- 
cuit to evaluate any function f{x) that can be evaluated classically). If the bit b is initially placed in a superposition 

-7=(|0) - |1)) , this circuit will invert the amplitudes precisely in the states for which f(x) = 1 , while leaving ampli- 

J2 

tudes in other states unchanged. 

4. Summary of Applications As mentioned in section 1, the search problem is the following: a function 
f(x), x = 0, 1 , . . . ( N - 1 ) , is given which is known to be non-zero at certain values of x ; the task is to find one such 
value. No structure is known for f(x) except for what is explicitly mentioned in the specific problems (4. 1 ) . . . (4.7) . 

n 

N is assumed to be a power of 2 , say N = 2 . There is a one-to-one correspondence between the N values of x 
and the respective states of an n-bit register. States corresponding to values of x for which f(x) is non-zero, are 
referred to as f-states. 

(4.1) Assume that f(x) = everywhere except for a single value of x. This is the standard problem of exhaustive 
search. 

(4.2) As in (4.1), there is a single point, t, where f(x) is non-zero. Some information about t is available in the fol- 
lowing form - another n bit word, r, is given which is known to differ from t in at most k out of the n bits. 

(4.3) There are multiple points (f-states) at which f(x) is non-zero, it is required to find any one of these. Some 
structure to the problem is specified in the following form. We are given a certain unitary transformation U & 



multiple s-states so that t/ for any s & any t are the same. The setting of this subsection is abstract, (4.4), 

(4.5) & (4.7) apply the framework of this subsection to concrete problems. 

(4.4) As in (4.3), there are multiple points at which f(x) is non-zero, it is required to find any one of these. How- 
ever, unlike (4.3), no further structure to the problem is given. 

(4.5) There is a single point t where f(x) is non-zero. Some information about t is available in the form of / n-bit 
strings, each of which differ from t in exactly k out of n bits. 

(4.6) f(x) = everywhere except at the unique point x = t , it is required to find t. Also, as in (4.3), we are given 
a unitary transformation U and multiple s-states. However, unlike (4.3), U f for various s-states & various t- 
states are not all identical. 

The analysis of section 2 extends the power of quantum search so that it can be used with an arbitrary unitary 
transform U , but only with a single s and single f-state. (4.3) extends it to multiple states, but in a restricted 
way. This derivation extends to multiple s-states. It is still not known how to handle multiple f-states that are 
not exactly symmetric. 

(4.7) This problem illustrates how the abstract techniques discussed earlier can be applied to solve an actual prob- 
lem. This problem was first discussed by Eddie Farhi & Sam Gutmann [Structure]. 

Two functions f(x, y) & g(x) are defined on the domain x = 0, 1, ...(N- 1) , y = 0, 1, ...(N- 1) . 
f(x, y) is zero everywhere except at the unique point (f p f 2 ) > g( x ) is non-zero at M values of x including 

x = ?j (here M « N ). The problem is to find fj & ? 2 . Classically this problem would take Cl(NM) steps. 

The algorithm of (4.1), without using the function g(x) , would take O(N) steps. The following analysis 

makes use of the general technique of (4.3) to develop an O(jNM) step algorithm. Several variants of this 
problem are also considered. 

4.0 The Approach The following general approach is made use of in each of the next 7 sub-sections - (4. 1 ) . . . (4.7) . 

There are N = 2" states, represented by n qubits, the task is to get the system into some target state(s) t at 
which f(x) is non-zero. A unitary transform U and the initial state s are selected and t/ is calculated. It then fol- 

71 -1 
lows by section 2 that by ^. , repetitions of the operation sequence -I g U I f U , followed by a single application 

I ts\ 

of U , the initial state s is transformed into the final state t. 

4.1 Exhaustive Search Assume that f{x) = everywhere except at a single point t. The object is to find t. 

As mentioned in the first paragraph of the introduction, there are several important problems in computer science 
for which there no solution is known, except exhaustive search. 



Solution For the W-H transform, described in section 3, U fs between any pair of states s & t is -~j= ■ Therefore we 

can start with any state s and the procedure of (4.0) gives us an algorithm requiring o( r— !— : j steps, i.e. O(Jn) 
steps. 

In case s be chosen to be the state, then the operation sequence Q = -I-WI W leads to the standard quantum 

1 

search algorithm based on the inversion about average interpretation [Gensrch] (note that W 1 = W ). To see this 

T ( T\ T 

write /- as I -2v-v- . Therefore for any vector x: -WI-Wx = -W \I-2v-v- \Wx = - x + 2Wv-v- Wx . It is 

00 } \ 0) 00 

T 

easily seen that Wv-v-Wi is another vector each of whose components is the same and equal to A where 
N-l 

A = — ^ x- (the average value of all components). Therefore the i th component of -WI-Wx is simply: 

1 = 

(-x. + 2 A) . This may be written as A + (A - x.) , i.e. each component is as much above (below) the average as it was 

initially below (above) the average, which is precisely the inversion about average. 

In case s be chosen to be a state different from , the dynamics is still very similar to the standard quantum 
search algorithm; however, the inversion about average interpretation no longer applies. 



4.2 Search when an item near the desired state is known: This problem is similar to (4.1), i.e. f(x) = except at 
the single point t. The difference from (4.1) is that some information about the solution, t, is available in the following 
form: another n bit word, r, is specified - 1 is known to differ from r in at most k of the n bits. 

Such a problem would occur in any situation when we had some prior information about the solution, this infor- 
mation could come either from prior knowledge or from a noisy data-transmission. 

Solution: The effect of the constraint is to reduce the size of the solution space. One way of making use of this con- 
straint, would be to map this to another problem and then exhaustively search the reduced space using (4.1). However, 
such a mapping would involve additional overhead. This section presents a different approach which carries over to 
more complicated situations as in (4.5). 

Instead of choosing U as the W-H transform, as in (4.1), in this section U is tailored to the problem under 
consideration. The starting state s is chosen to be the specified word r. The operation U consists of the following uni- 



tary transformation 
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applied to each of the n qubits. Calculating U fs , it follows that 



n — k k 

\U f \ = fl--l f-1 and log It/, I = ^log- — --^log^—^. The technique described in (4.0) can now be 

\ ts \ \ nj \nj I "I 2 n 2 k 

used - as in (4.1), this consists of repeating the sequence of operations -/ s UI f U , ^(p" times, followed by a sin- 
gle application of the operation U (note that, as in (4.1), U 1 = U ). 

( ll\ n ' 

The size of the space being searched in this problem is approximately LI which is equal to - — . Using 

( Yl\ k 

Stirling's approximation: logn! = nlogn - n , it follows that log , ~ nlog - felog , comparing this to the 

\kj n-k n-k 

number of steps required by the algorithm, we find that the number of steps in this algorithm, as in (4.1), varies as the 
square-root of the size of the solution space being searched. 



4.3 Multiple s & t states with the same t/ : f(x) is non-zero at P values of x, i.e. there are p f-states. Some 

structure of the problem is specified in the following form. Assume that we have at our disposal a unitary transform 
U and a s-states such that U {s between any f-state and any s-state is the same. The object is to find one of the t- 

states. This is accomplished by transforming the system into a superposition so that there is an equal amplitude in 
each of the f-states and zero amplitude elsewhere. After this, a measurement is made that projects the system into one 
of its basis states, this gives a f-state. 

The problem considered in this subsection is abstract in the sense U is an arbitrary unitary transformation. (4.4), 
(4.5) and (4.7) apply this to concrete problems. 

Solution: The approach is similar to the exhaustive search problem of (4.1). However, the analysis of section 2 has to 
be redone with the following three changes: 

a- 1 

1 

(a) The starting state instead of being v , is the superposition — } v„ - the amplitude in all s states is equal 

4 ' Ja ^ s ■• " 

a = 

to -|= , and zero everywhere else. Assuming a to be a power of 2 ( a = 2° ), such a superposition can be easily 

created by the following procedure. Start with an a bit system with all bits in the state. Do a W-H transform 

on the a bit system and then carry out a mapping from the 2° states to the s-states. 

(b) The operations I ^ & / f invert the amplitudes in all s-states & all f-states, respectively. 



(c) It can then be shown by an analysis similar to section 2, that after O 



1 



operations of -I \J '/ { U 



followed by a single application of U , the system reaches a superposition so that the amplitude is equal in all 



(3 f-states and is zero everywhere else. Note that the number of operations is smaller by a factor Jaft as com- 
pared to the situation when there were single s-states & single f-states (as in (4.1)). 



4.4 Problem f(x) is non-zero at (3 values of x , equivalently there are P f-states - the task is to find one of these. 

This is the problem of exhaustive search when there are multiple (P) solutions. A classical search would take 

an average of O^jj^j steps to find a solution. This section presents an J^j step quantum mechanical algorithm. 
Solution By the definition of the W-H transform in section 3, W - for any f is the same. Therefore if we choose s as 

the state, then it follows by (4.3) that after repetitions of -I s WI ( W followed by a single application of 

W, the system reaches a superposition such that the amplitude is equal in all the f states and zero everywhere else. 
Note the following three points regarding this scheme: 

As in (4.3), the operation inverts the phase for all P f-states. 

The above implementation requires P to be known in advance. 

The search time is VP faster than the exhaustive search algorithm of (4.1). 

It is necessary to choose s as the state, this is different from (4.1) where s could be arbitrary. 

4.5 Problem f(x) = except at the single point f. Some information about f is available in the form of a rc-bit 
strings, each of which differs from f in exactly k bits. 

This is in some sense the dual of (4.4). In that case there were multiple f-states but a single s-state, while in this 
problem there are multiple s-states and a single f-state. This kind of problem could occur in extracting a signal out of 
multiple noisy transmissions. 

Solution Let the a specified states be the s-states. Initialize the system to a superposition of these states by the pro- 
cess described in (4.3)(a). After this, apply the unitary transform U which applies the following unitary operation 
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n Nn 



n 



n-k k 

to each qubit. As in (4.2), It/. I = I 1 - -1 f-~\ and log l/J. I = ^log- — - - ^log^-^ for 

I "I V nj ynj I "I 2 n 2 k 



all i-states. Also, since each of the s-states differ from f in exactly the same number of bits implies that U fs has the 

same sign for all s-states. The framework of (4.3) can now be used - this yields an algorithm that is times faster 
than that of (4.2). 



4.6 Multiple s-states & a single <-state such that t/ between various s-states & t is not identical f(x) = 

everywhere, except at the single point t. As in (4.3), some structure to the problem is specified via U . Assume that we 
have at our disposal a unitary transform U and various s-states are specified. However, unlike (4.3), 17 between 

various s and various f-states are not exactly equal. This case is qualitatively different from all of those considered so 
far because the analysis of section 2 or the modified analysis of (4.3), does not directly apply. 

This extends the noisy data-transmission problem of (4.5), to the case where there are multiple states specified 
that are close to the solution state but differ from it in varying number of bits ((4.5) required each of the given states 
to differ from the solution in exactly the same number of bits). 

Solution Assume the number of ^-states to be a , further assume that a is a power of 2, i.e. a = 2° . Consider V to 

-1 -1 -1 -1 

be a unitary matrix which is a product of 3 simpler unitary matrices, i.e. V = V j V 2 ^ an d ^ = U V2 . 
Vj is a W-H transformation on a bits and a states where a = 2° , 

V 2 maps the a = 2 a states generated by V j onto the respective s-states in jV-dimensional state space, 
U is the available unitary transform on the N states. 

Let the initial state be the state. As a result of V j V 2 , the amplitude in each of the states: Sq, s^...s 



becomes — ; after U , the amplitude in the f-state is — ) J7 . By (4.0), it follows that after ; repeti- 

Ja ' Toe ^ a a - 1 

a = 



a = 



ts, 



tions of Q = -I- V I V followed by a single application of V , the amplitude in the target state becomes O(l). 
' 

Many of the problems in the previous subsections can be seen to be particular cases of this. For example in case 

U are all equal, say to u , then the number of iterations becomes: — ^— . The algorithm and bound of (4.5) imme- 
a \u\Ja 

diately follow from this. 

4.7 Two dimensional search Two functions f(x,y) & g(x) are defined on the domain x = 0, 1, ...(N - 1) , 
y = 0, 1, ...(N - 1) . f(x, y) is zero everywhere except at the unique point (fp ? 2 ) , g(x) is non-zero at M values 

of x including x = t j (here M « N ). The problem is to find t j & ? 2 . 

Classically this problem would take Q.(NM) steps. The algorithm of (4.1), without using the function g(x) 

would take O(N) steps. This section presents an O(jNM) step algorithm. For a different analysis, along with a 



proof that the algorithm of this section is within a constant factor of the fastest possible algorithm, see [Structure]. 

Several variations of this problem are also briefly considered, these demonstrate how the techniques discussed 
in this paper can be applied to real problems. 

Solution First consider the function g(x) , x = 0, 1, ... (N - 1 ) . By executing the algorithm of (4.4), with the f-states 
as the non-zero values of g(x) , it is possible for the system to reach a superposition such that at each point at which 

g(x) is non-zero, the amplitude is —Lr and the amplitude is zero everywhere else. This is accomplished by a 

Jm 

sequence of J^j unitary transformations (4.4) - denote this composite unitary operation by t/j . Next keep the 

value of x the same and carry out the algorithm of (4. 1) on the N values of y with the f-state corresponding to the non- 
zero value of the function f(x, y) . This consists of a sequence of 0(J~N) elementary unitary transformations, 
denote this composite unitary operation by t/ 2 ■ 11 follows from (4.1), that as a result of this operation sequence, in 

case the particle is at x = 1 1 , it is also at y = ? 2 . The amplitude of the system being in the desired state is therefore 
— != . By means of the unitary transformation U = U, f/ 9 , the system starting from a certain initial state reaches the 

Jm 1 1 

desired state with an amplitude of — != . 

Jm 

Since U is a sequence of elementary unitary operations, it follows that U 1 is a sequence of the adjoints of 
the same operations in the opposite order and can hence be synthesized. By applying the procedure of section 2, it fol- 
lows that by repeating the sequence of operations Q = -I U I ( U , O(jM) times, the system reaches the target state 
with certainty. 

The total number of steps is given by the number of repetitions of Q (i.e. 0(jM) ) times the number of steps 

required for each repetition, (i.e. {o{^ J^j + 0( jN)j ) which gives 0(jM) x ^ J^j + 0( jN)j = O( jNM) . 

The above analysis easily extends to more general cases. For example consider the case, where the search- 
space is rectangular instead of square, i.e. the number of possible values for x is JVj and for y is N 2 , The number of 

steps, instead of being O(jNM) , now becomes 0(jN~^) + 0( Jn^M) . Alternatively consider the case where there 
is an r) dimensional space with the same number of points (N) in each dimension. Instead of just f(x, y) & g(x) , 
there are now r) functions: f(x^, x^ ■■■x^) , g^x^x^ ■■■ x r ^ _ i ) > ^2^ x l' x 2' "" x n ~t)' S^_\i x \) w i tn analo- 
gous definitions. It follows by a similar approach that the number of steps is now 0(jNM jMj.-.M^ _ j) . 

Another variation is when the function f(x, y) is non-zero at multiple points, say (3 points, and it is required 



to find one of these. In case each of the P points has a different value of x , the framework of (4.3) applies. Consider- 



ing the unitary transformation U as U^U^, leads to an algorithm that requires ^"j^J steps. 

It is not clear how to derive an algorithm when some of the P points have the same value of x . 

5. Conclusion & Further Work [Search] demonstrated how to make use of quantum mechanical properties to 
develop exhaustive search kinds of algorithms, i.e. algorithms for problems that lacked any structure. [Search] used 
subtle properties of a particular quantum operation called the W-H transform. Subsequently [Gensrch] extended this 
so that other quantum operations could be used instead of the W-H transform. 

Most interesting problems in computer science are concerned with the structure of problems and how to 
develop algorithms to take advantage of this structure. [Hogg] has previously suggested heuristic quantum mechani- 
cal algorithms for structured problems. This paper has given several examples of structured problems and how 
search-type algorithms can be extended to solve these ((4.2) through (4.7)) - quantitative closed form bounds were 
derived for the running time of these algorithms. The extensions have shown how to deal with the situation where 
there are multiple ^-states (initial states) and a single f-state (target state) ((4.5) and (4.6)). Also, it is possible to deal 
with multiple f-states provided they are exactly symmetric ((4.3) and (4.4)). The next step would be to obtain a gen- 
eral algorithm with multiple f-states. 
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